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Abstract 



We analyze the coherence properties of polarized neutrons, after they have interacted with a magnetic field or a phase shifter 
undergoing different kinds of statistical fluctuations. We endeavor to probe the degree of disorder of the distribution of the 
phase shifts by means of the loss of quantum mechanical coherence of the neutron. We find that the notion of entropy of 
the shifts and that of decoherence of the neutron do not necessarily agree. In some cases the neutron wave function is more 
coherent, even though it has interacted with a more disordered medium. 



I. INTRODUCTION 



The notion of dccoherence has attracted increasing attention in the literature of the last few years pJ2] . The loss 
of quantum mechanical coherence undergone by a quantum system, as a consequence of its interaction with a given 
environment, can be discussed in relation to many different physical phenomena and has deepened our comprehension 
of fundamental issues, disclosing unexpected applications as well as innovative technology. 

Neutron physics (neutron optics in particular) has played an important role in this context, both on theoretical 
and experimental grounds. Non-classical states are readily obtained, for instance by splitting and then superposing 
wavepackets in an interferometer || or different spin states in a magnetic field j||| , and are of great significance in the 
investigation of fundamental quantum mechanical properties. The aim of this paper is to investigate the coherence 
features of neutron wave packets, by making use of the Wigner function |6|, in analogy with concepts and techniques 
that are routinely used in quantum optics JM] . The studies of the last few years have shown that non-classical states 
are fragile against statistical fluctuations [g[P] : the analysis of situations where these states display robustness during 
the interaction with noisy environments is therefore of great practical interest. 

The main motivation of this work is to use the coherence properties of the wave function as a "probe" to check 
the degree of disorder of an environment. A similar idea was first proposed, as far as we know, in the context of 
quantum chaos and Feynman integrals Jl0| . One might naively expect that a neutron ensemble suffers a greater 
loss of quantum coherence by interacting with an increasingly disordered environment: intuitively, a more disordered 
environment should provoke more randomization of the phase of the wave function, which in turn implies more 
quantum decoherence. As we shall see, this is not always true: some of the results to be discussed below are rather 
counterintuitive and at variance with naive expectation. In some cases the neutron wave function is more coherent, 
even though it has interacted with a more disordered medium. This statement can be given a precise quantitative 
meaning in terms of the entropy of the medium and of a "decoherence parameter" that will be defined for the neutron 
density matrix. 



II. PRELIMINARIES 



The Wigner quasidistribution function || can be defined in terms of the density matrix p as 

W(x, k) = ^j<% e~ ik Hx + £/2\p\x - £/2), 



(1) 



where x and p — hk are the position and momentum of the particle. One easily checks that the Wigner function is 
normalized to unity and its marginals represent the position and momentum distributions 



Trp= dxdkW(x,k) = 1, 
P(x) = (x\p\x) = I dk W(x,k), 
P{k) = (k\p\k) = J dxW{x,k). 



Notice that 



dx dk W(x, k) 
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In this paper we shall consider a one-dimensional system (the extension to 3 dimensions is straightforward) and assume 
that the wave function is well approximated by a Gaussian 



ip(x) = (x\ip) = 
0(fc) = (A#> = 
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where ip{x) and (f>(k) are the wave functions in the position and momentum representation, respectively, 8 is the 
spatial spread of the wave packet, SkS = 1/2, Xo is the initial average position of the particle and pa = ftfco its average 
momentum. The two functions above are both normalized to one: normalization will play an important role in our 
analysis and will never be neglected. The Wigner function for the state (||)-(0) is readily calculated 



W(x, k) = — exp 

IT 



(x-x ) 5 



2<5 2 



cxp[-2S 2 (k-k ) 2 ] • 



(8) 



In this paper we will focus on two physical situations. In the first one, a polarized neutron acquires a phase shift 
A, either by going through a phase shifter or by crossing a magnetic field parallel to its spin. In the second one, a 
polarized neutron is divided in two states, either in an interferometer or by crossing a magnetic field perpendicular 
to its spin. The latter situation is physically most interesting, for it yields non-classical states, whose coherence 
properties are of great interest. 

A. Single Gaussian 

If a Gaussian wave packet undergoes a phase shift A, the resulting Wigner function reads 



1 



W(x, k, A) = — exp 



(x - x + A) 2 



25 2 



exp[-2(5 2 (fc-fco) 2 ] 



(9) 



Physically, this is achieved either by placing a phase shifter in the neutron path, or by injecting a polarized neutron 
in a constant magnetic field parallel to its spin. In both cases, the total energy of the neutron is conserved. In the 
latter case, if the field has intensity B and is contained in a region of length L, the neutron kinetic energy in the field 
changes by AE = —\/i\B, where |/Lt| is the neutron magnetic moment. This entails a change in average momentum 
Afc = m^,B/h 2 ko and a phase shift proportional to A = LAk/ka. When it leaves the field, the neutron acquires again 
the initial kinetic energy. 

B. Double Gaussian 

Consider now a neutron wave packet that is split and then recombined in an interferometer, with a phase shifter 
placed in one of the two routes. The Wigner function in the ordinary channel (transmitted component) is readily 
computed: 



W mt (x,k,A) 
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cos(fcA) 



(10) 



Notice that, for A ^ 0, it is not normalized to unity (some neutrons end up in the extraordinary channel — reflected 
component) and that for A = (no phase shifter) one recovers (S). 

A similar result is obtained when a polarized (say, +y) neutron crosses a magnetic field aligned along an orthogonal 
direction (say, +z). The total neutron energy is conserved, but due to Zeeman splitting the two spin states in the 
direction of the B field have different kinetic energies and travel with different speeds. This is a situation typically 
encountered in the so-called longitudinal Stern-Gerlach effect 0] and in neutron spin-echo experiments m (except 
that we are not considering the second half of the evolution, with an opposite B field that recombines the two spin 
states) . An experimental realization of this situation was investigated very recently pdj] . If the initial wave function 



l*> = IVO 



\lp) <g> 
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—I 

V2 l 



(11) 



where |±) Q (a = x,y,z) represents spin up/down in direction a, the final state in the position representation, after 
crossing the S-field, reads 



(af |*>- ^^ (« + f ) |+). + -^ 



ip I a; 



(12) 



If only the +y-spin component is observed ("post selection" of the initial spin component jl2[]) the probability ampli- 
tude is 



v (+,x\*) = - 

and the Wigncr function is readily computed as 

1 



i/j [x + — ] +ip [x 



(13) 



W magn (x,/c,A) 
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cos(fcA) 



(14) 



This result is slightly different from (|10[), because in this case both spin components undergo a phase shift (±A/2). 
Once again, for A = (no magnetic field) one reobtains (pi). 

We stress that in both cases the neutron wave packet has a natural spread 6t = \/ S 1 + (Ht/2mS) 2 (due to its free 
evolution for a time t ~ mL/hko); however, this additional effect will be neglected, because, as proved in Appendix 
A, it is not relevant for the loss of quantum coherence. 



III. FLUCTUATING PHASE SHIFT 



The previous analysis refers to a rather idealized case, in which every neutron in the beam acquires a constant 
phase shift. This is clearly not a realistic situation, for it does not take into account the statistical fluctuations of the 
B field or of the shifter in the transverse section of the beam. If, for any reason, the phase shift A fluctuates, the 
neutron beam will partially loose its quantum coherence and the Wigner function will be affected accordingly. We 
shall consider the case of "slow" fluctuations, in the sense that each neutron crosses an approximately static B field 
(or a phase shifter of uniform length L) , but the intensity of the field (or the length of the shifter) varies for different 
neutrons in the beam (different "events" ) . We will suppose that every neutron undergoes a shift A that is statistically 
distributed according to a distribution law w(A). The collective "degree of disorder" of the shifts A can be given a 
quantitative meaning in terms of the entropy 



S = - f dAw(A) log(w(A)). 
On the other hand, the average Wigner function reads 



W m (x,k) 



f dAw(A) W(x,k,A) 



(15) 



(16) 



and represents a partially mixed state. The coherence properties of the neutron ensemble can be analyzed in terms 
of a decoherence parameter ]l3[ | 



Trp 2 _ 2irfdx dk W m (x, kf 
(Trp) 2 (JdxdkW m (x,k)) 2 



(17) 



This quantity measures the degree of "purity" of a quantum state: it is maximum when the state is maximally mixed 
(Trp 2 <Trp) and vanishes when the state is pure (Trp 2 =Trp): in the former case the fluctuations of A are large and 
the quantum mechanical coherence is completely lost, while in the latter case A does not fluctuate and the quantum 
mechanical coherence is perfectly preserved. The parameter (|17| ) was introduced within the framework of the so-called 
"many Hilbert space" theory of quantum measurements fipfand yields a quantitative estimate of decoherence. The 
related quantity Trp— Trp 2 (that might be called "idempotency defect") was first considered by Watanabe |1J] many 
years ago. A measure of information for a quantum system has been recently introduced, which is related to e and is 
more suitable than the Shannon entropy JL5| . 

One might naively think that the two quantities S and e should at least qualitatively agree: in other words, the loss 
of quantum mechanical coherence should be larger when the neutron beam interacts with fluctuating shifts of larger 
entropy. Such a naive expectation turns out to be incorrect. Our purpose is to investigate this problem. To this end, 
it is useful to consider some particular cases. 



A. Gaussian noise 



We first assume that the shifts A fluctuate around their average Aq according to a Gaussian law: 



iu(A) = 



1 



V27R7 2 



exp 



(A-A„) s 



2a 2 



(18) 



where a is the standard deviation. The ratio ct/Ao is simply equal to the ratio SB /Bo (or SL/Lq), SB (SL) being the 
standard deviation of the fluctuating magnetic field (length of phase shifter) and Bq (Lq) its average. The entropy of 
( |l8| ) is readily computed from ( |lq ) 



and is obviously an increasing function of a. 



S = -log(27rea 2 ) 



1. Single Gaussian 



(19) 



Consider now a neutron described by a Gaussian wave packet. If the phase shift A fluctuates according to (jig), 
the average Wigner function is readily computed by (|6|), (^ and (jig), 



S 2 



W m (x,k) 



and its marginals (||)-(§) are easily evaluated 



P(x) = 



exp[-25 2 (k - k ) 2 } exp 



(x - x + A ) 2 



2{8 2 + a 2 ) 



^/2n(5 2 + a 2 ) 

^s 2 



exp 



(X -XQ + Ap) ; 

2(S 2 + a 2 ) 



P(k) = \ exp [-2S 2 (k - fc ) 2 ] . 



(20) 



(21) 



(22) 



Notice that the momentum distribution (E2) is unaltered and identical to |</>(A:)| 2 in (ffl): obviously, the energy of each 
neutron does not change. Observe on the other hand the additional spread in position S' = \/5 2 + a 2 (Figure |l|) and 
notice that the Wigner function and its marginals are always normalized to one. 
The decoherence parameter (uj\) can be analytically evaluated 



e= 1- 



S 2 



S 2 + a 2 



(23) 



and is a monotonic function of a for every value of 5. This behavior is in qualitative agreement with that of the 
entropy (fl9) . As expected, a more entropic distribution of phase shifts entails a greater loss of quantum mechanical 
coherence for the neutron ensemble. The behavior of e vs S and a is shown in Figure 0(a). 



2. Double Gaussian in an interferometer 



Consider now the double Gaussian state (|10|), obtained when a neutron beam crosses an interferometer. The average 
Wigner function ( |l6| ) reads 



W^\x,k) 
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Ait 
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+ exp 



x 
'2S 2 



(24) 



where we set xq — for simplicity. Its marginals 
momentum probability distribution reads 



2{S 2 + ^) )Y 
can both be computed analytically; in particular, the 



P(fc) = \/^exp[-2<5 2 (fc-fc ) 2 ] 



1 + cxp - 



cos(fcAo) 



(25) 



As one can see from Figure ia) , interference is exponentially suppressed at high values of k and the oscillating part 
of the Wigner function is bent towards the negative s-axis. This is due to the ^-dependence of the cosine term in 
(p4|) that entails different frequencies for different values of x. The decoherence parameter (p7|) reads 




exp 



2<5 2 
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S 2 



-k Ao 



(27) 



represents the probability of detecting a neutron in the ordinary channel. The explicit expression Q2q ) of the deco- 
herence parameter is involved and difficult to understand. Therefore, e is shown in Figure 0(b) as a function of 5 and 
a for fixed values of fco and Ao: somewhat surprisingly, for some values of 5, even though the noise a increases, the 
decoherence e decreases. 

Observe also that e never reaches unity: e < 3/4. This is due to the fact that one of the two Gaussians does not 
undergo any fluctuations (there is a fluctuating phase shifter in only one of the two routes of the interferometer): 
therefore a part of the Wigner function is not affected by noise, as one can see in Figure 0(a). We shall comment 
again on the peculiar features of e in a while. 



3. Double Gaussian in a magnetic field 

If we consider a polarized neutron beam interacting with a B field perpendicular to its spin, Eq. ( |l4| ) yields 

exp[-2<5 2 (fc-fc )2i 



WC agn (zAA) 
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S 2 + ^ 



exp 
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2(J2 + 21) 




2 exp — — -p? 



(28) 



This Wigner function has the same k- marginal ( |25| ) as the previous one (although the rr-marginals are different) . Also 
in this case, one observes a strong suppression of interference at large values of momentum P.O.f3[, but without 
cc-dependence in the cosine. See Figure |](b). In this case, the decoherence parameter ( |f7| ) reads 
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(29) 



Again, the explicit expression of the decoherence parameter is complicated and depends on several physical parameters; 
it is therefore convenient to concentrate on a particular case. An experimental realization of a fluc tuati ng shift 



(according to a given statistical law) is easier with the magnetic field arrangement discussed in Section II B . Let us 
therefore consider the experiment |ll[ |, in which a polarized (+y) neutron enters a magnetic field, perpendicular to 
its spin, of intensity Bq = 0.28mT, confined in a region of length L = 57cm. The average neutron wavenumber is 
fco = 1.7 • 10 10 m _1 and its coherence length (defined by a chopper) is 6 = 1.1 • 10~ 10 m. By travelling in the magnetic 
field, the two neutron spin states are separated by a distance Ao = 2miiBoL/ti 2 kQ — 16.1 • 10~ 10 m, one order of 
magnitude larger than <5. The behavior of e in (|29J) is shown in Figure 2(c) for these experimental values: observe 
that for S > 3A e is not a monotonic function of a: in other words, for some values of the parameters, even though 
the noise a increases, the decoherence e decreases. This is at variance with the behavior of the entropy ([19]) and with 
what one might naively expect. We conclude that, in general, both for a double Gaussian in an interferometer and in 
a magnetic field, the behavior of e does not agree with that of the entropy. 

B. Sinusoidal fluctuations with increasingly less rational frequencies 

In order to shed some more light on the results of the previous subsection, let us consider a different example, that 
is more convenient for an experimental perspective. Suppose that the phase shift changes according to the law 

A(£) = A + Ai[sin(fit) + sinfofi*)]. (30) 

where t is time, 51 a frequency much smaller than Vq/L, the inverse time of flight of the neutron in the shifter's region, 
Ao the mean phase shift, Ai(<C Ao) the "fluctuation" width (see below) and rj (0 < r? < 1) a real number. For the 
neutron ensemble (the beam) the shifts will be distributed according to law 

«;(A)= fdtf(t) 6(A-A(t)), (31) 

where f(t) is the probability density function of the stochastic variable t. In our case, f(t) = 1/T in (0,T), where 
T(3> SI -1 ) is a sufficiently large time interval. In such a case, by making use of (pl|), the Wigner function can be 
expressed as an ergodic average 

W m {x, k)= I dA w(A) W(x, k,A) = - J dt W(x, k, A(t)). (T large) (32) 

We stress that A is treated like a random variable although, strictly speaking, the underlying process is deterministic. 
However, this is not a conceptual difficulty: in practice, one just treats the neutron ensemble in an experimental run 
without looking at the correlations among different neutrons. The same effects on the neutron ensemble would be 
obtained by first generating a random variable t, uniformly distributed in (0,T), then constructing the additional 
random variable A according to ( J30[ ) and finally accumulating all neutrons in the experimental run. In this way 
different neutrons are uncorrelated. The distribution law of the shifts (|3(j) can be obtained by means of a i3-field 

B{t) = B + Bi [sin(fit) + sinfoOt)] . (33) 

Like in the previous section, we assume that B is a slowly varying function of time, so that each neutron experiences 



a static field during its interaction. Observe that the scheme proposed in (J30|)- ( 32 ) is not difficult to realize experi- 
mentally. On the other hand, it would be complicated to obtain the same distribution of shifts with a phase shifter 
placed in one of the two routes of an interferometer. 

We will study the coherence properties of the neutron beam when it crosses a magnetic field made up of two 
"increasingly less rational" frequencies, by choosing 

— fj (34) 



J 3+1 

where fj are the Fibonacci numbers 

fi+i = fj + f j-i (/o = /i = l). (35) 



This particular choice is motivated by the (naive) expectation that an oscillating magnetic field (33) composed of 



mutually less rational frequencies should provoke more decoherence on the neutron ensemble. Once again, this 



3- 



expectation will turn out to be incorrect. The ratios ( pif ) tend to the golden mean (the "most irrational" number 
J) as j increases 

In general, one cannot obtain an analytic expression for the probability density function (Bh); however, an accurate 
numerical evaluation of iu(A) is possible: for every finite value of j, r.j is a rational number, so that one can integrate 
( pl| ) over the interval T = Jj + i2it /Vt. In Figure |4| we show the results of our numerical analysis. The distribution 
function w(A) has a finite number of (integrable) divergences in its interval of definition; as the order in the Fibonacci 
sequence becomes higher, the number of divergences in the interval grows. In the j — ► oo limit, i.e. for the golden 
mean r^ = (V5 — l)/2, it is possible to apply the theorem on averages for the ergodic motion on a torus ||17| and 
find an analytical expression of w(A) in terms of an elliptic integral of first kind (see Appendix B). The resulting 
distribution is a smooth function with only one (integrable) divergence in A = and is plotted in Fig. |f|(f). 

By applying the same technique utilized for the numerical evaluation of w, the entropy is computed according to 
the formula 

S = -JdA w(A) \og(w(A)) =~J dt logHA(t))], (37) 

which is easily obtained by Eqs. ( |i"5| ) and (3l|) (using the value T — /j + i27r/f2 for the numerical evaluation). 

The decoherence parameter is computed from Eq. (17), first with the Wigner function (|l) (single Gaussian) and 



then with the Wigner function (pTJ) (double Gaussian in a magnetic field): in both formulas, we used Eq. ( p2| ) and 
set Ao = 16.1 • 10 _10 m, Ai = 2 • 10~ 10 m and the same numerical values of the previous subsection for fco and 5 [ pT[ . 
Our results are summarized in Tabic 1 and Figure ]a. 

We notice that, although, for j = 1, ...,5, S is a monotonically increasing function of the Fibonacci number in the 
sequence, e reaches a maximum for Tj — 3/5 (i.e., j = 3). It is remarkable that the maximum is obtained for the same 
Fibonacci ratio in both cases (single and double Gaussian). Once again, the behavior of entropy and decoherence are 
qualitatively different. Figure 5| should be compared to Figure |2j: it is worth noting that in the case analyzed in this 

, the behavior of entropy and decoherence do not agree even when the neutron state is 



section, unlike in Section III A 



a single Gaussian (namely, a "classical" state). 

IV. CONCLUSIONS 

Decoherence is a very useful concept, that has recently been widely investigated and has turned out to be very 
prolific. It is intuitively related to the loss of "purity" of a quantum mechanical state and can be given a quantitative 
definition, as in (fL7j). However, we have seen that the very notion of decoherence is delicate: in particular, it is not 
correct to think that a quantum system, by interacting with an increasingly "disordered" environment, will suffer an 
increasing loss of quantum coherence. Our analysis has been performed by assuming that each neutron, during an 
experimental run, interacts with a constant magnetic field: the neutron beam, on the average, undegoes decoherence. 
This "quasistatic" ap proxim ation is only a working hypothesis and will be relaxed in future work. As emphasized at 



the beginning of Sec. IIIB, it is easy to achieve experimentally in the limit Q -C vq/L, where Q is a characteristic 
frequency of the fluctuation and (vo/L)^ 1 the time of flight of the neutron in the phase shifter or in the magnetic 
field. This approximation also simplifies (both conceptually and technically) our theoretical analysis, without however 
having a substantial influence on our general conclusions. It is worth stressing, in this respect, that the decoherence 
parameter, defined in fll7|), depends on the interaction and not on the free Hamiltonian, at least in the physical 
situations investigated here (see Appendix A). 

Our quantitative definition of decoherence depends, as it should, on the very characteristics of the experimental 
setup: the decoherence parameter is defined in terms of the average Wigner function (or equivalcntly the density 
matrix) of the neutron ensemble, after the interaction with the apparatus. An experimental check of the features 
of the Wigner functions discussed in this paper would require its tomographic observation. Similar techniques are 
commonly applied in quantum optics pM and wou ld be available in neutron optics as well, in particular for the 



experimental arrangement discussed in Section III E . However, we think that a better comprehension of the effects 
analyzed in this paper could probably be achieved by studying the marginals (or possibly some other tomographic 
projection) of the Wigner function and the visibility of the interference pattern. Additional work is in progress in 
this direction. From an experimental perspective, an analysis of decoherence effects along the guidelines discussed 
here would be challenging: although the concepts of decoherence and entropy are intuitively related, they display 
some interesting differences. If properly understood, those situations in which a larger noise yields a more coherent 
quantum ensemble might lead to unexpected applications. 
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APPENDIX A 



We prove that the decoherence parameter wj\), under rather general conditions, does not depend on the free 
evolution of a quantum system. Let the Hamiltonian of a quantum system be 



H = H a + H 1 (a), 



(A.l) 



where Hq and Hi are the free and interaction Hamiltonians, respectively, and a is a c- number (that can fluctuate 
according to a given statistical law). We assume that 



[H ,Hi{a)} = iC, with [H , C] = [Hi(a), C] = 0, 
where C is in general an operator (independent of a), so that (h = 1) 



e -it(H +Hi[a)) _ e it 2 C/2 e -itHo e -itH 1 (a) 



Consider now the density matrix at time t, 



Pa{t) 



_ p -it{H +H l ( a )) 



Po e 



it(J?o+fl"i(a)) 



where po is the initial density matrix. From Eq. (A.3) 



p a (t) = e lt2c ' 2 e~ ltHa e - ltH ^ p e ltH ^ e *n» e -i?C/2 
and the average over a yields 

where w(a) is the distribution function and the bar denotes average. Therefore 

Tr 



p(t) 



= Tr 



Pint(t) 



where p- ln t is the density matrix in the following interaction picture: 

p int (t) = e *Ho^t 2 C/2 p(f) e ^tH 0+l t*C/2 = e -UH x { a ) pQ e U Hl (a) 

This proves that the trace of the average density matrix does not depend on the free evolution. The result ( |A.7| ) can 
be generalized to any function of the average density matrix 



(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 

(A.7) 
(A.8) 



Tr 



f(p(tj) 



Tr 



/ Pint(t) 



This shows that the decoherence parameter defined in fll7|) does not depend on the free evolution: 



£int; 



(A.9) 



(A.10) 

as claimed at the end of Sec. ||. This result can be applied to the case studied in Sec. [II, where Hq = ■&—, Hi = —fiB 
and the parameter a is the intensity of the magnetic field B (whose direction is supposed constant). Notice that we 
are considering wave packets that interact with a constant and homogeneous field B from the initial time t = to the 
final time t ~ mL/hko, so that condition ( |A.2| ) is fulfilled. The case of a neutron wave packet in an interferometer is 
analogous, if we assume that the phase shifter simply yields a phase (optical potential approximation). 



APPENDIX B 



We compute here the distribution function (|3l| ) when f(t) = \jT and A(t) changes according to (30), in the j = oo 
limit (p6|). The function w(A) has a finite number of (integrable) divergences in its interval of definition. Notice 
that, as the order in the Fibonacci sequence becomes higher, the number of divergences in the interval grows and the 
numerical evaluation of the distribution function becomes more difficult. Let us introduce the two-component vector 

<Pj(t) — (<pi, if2) — Wjt, where Uj — (CI, rjfl). (B.l) 

The vector cp, performs a (quasi)periodic motion on the two-dimensional torus T 2 . In particular, for every finite 
value of j the frequencies are dependent, i.e. o^/wi = rj G Z, and the orbits are closed. For larger values of j the 
number of windings in a period increases and the length of the periodic orbit becomes larger. In the j — > 00 limit, 
the two frequencies become independent and the resulting motion on the 2-torus becomes ergodic: the trajectory is 
everywhere dense and uniformly distributed on T 2 . In this case, according to the theorem on averages [O, the time 
average of every integrable function f(<p) (where tp = ip^) coincides with its space average, i.e. 

1 pT -, i-2tt i-2-ix 



lim - / dt f(<p(t)) = —^ I dipt / d<p 2 f(<p). (B.2) 

T^oo T J (2ny J J 



Applying Eq. (|B.2|) to the function 



/( V ) = 5(A-A 1 [sin^ 1 +sinH) (B.3) 



we obtain 



1 f T 

■;(A) = lim - / dt S(A - Ai[sin(Ot) + sin(r 00 fii])) 

= , n x 7 / dipx / dtp 2 S(A - Ai[sin<pi +smip 2 ]) 
(2'kY J q Jo 



= aT / ds Ps{s) Ps ^ A / Al _ s )> ( B - 4 ) 



where Ps is the sine distribution 



1 



2 ? r 



Ps(s) = tt dcp6(s- sin <p) = (B.5) 

2vr Jo 7TV1 - s 2 



After some algebraic manipulation one finds 



2 „ / 1 I ( A N 2 



2Ai 



tt,(A) = ^Ai J!, l arcBin ^TM , v 1 "l2ArJ ! (R6) 



where F(0, 7) is the elliptic integral of first kind |Hj] 







1 



F(/3, 7 )=/ da - (B.7) 

Jo y 1 — 7^ sin a 



The limiting distribution function (BJ3) is plotted in Fig. ^(f). 

Observe that in Fig. ftfa-f) the number of divergences increases so quickly that, in the j = 00 limit (golden mean), 
w(A) becomes a smooth function with only one (integrable) divergence in A = (indeed w(A) ~ log(8Ai/|A|)/7r 2 Ai 
for A — ► 0) . In this sense Berry et al. coined the epigram "stocasticity is the ubiquity of catastrophe" pQ| . (Inciden- 
tally, notice the similarity of Fig. Q with Fig. 12 of |20[| .) 
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Fluctuating 
phase shift 
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FIG. 1. If the phase shift fluctuates, each wave packet acquires a different shift. This is pictorially represented in the figure, 
where different outgoing wave packets are displayed, each associated with a single neutron ("event"). The average Wigner 
function is given by Eq. (E0J). 
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FIG. 2. Decoherence parameter vs coherence length of the wave packet 5(A) and standard deviation of the fluctuation cr(A). 
(a) Gaussian wave packet, (b) Double Gaussian in an interferometer, with Aq = 16.1 A. (c) Double Gaussian in a magnetic 
field, with Ao = 16. 1A. In all cases fco = 1.7A -1 . Observe that in case (a) the decoherence parameter is a monotonic function 
of a for every value of S, while this is not true for cases (b) and (c). Notice also that in case (b) the decoherence parameter 
never reaches unity (e < 3/4): this is due to the fact that only one Gaussian (in one branch of the interferometer) undergoes 
statistical fluctuations [see Figure 0(a)]. 




(a) 



(b) 
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FIG. 3. Wigner functions for different values of the standard deviation a in (tlq). (a) Double Gaussian in an interferometer 
(E4j). (b) Double Gaussian in a magnetic field (|28[), From top to bottom, a — 0, 0.6, 1.2, 1.8A. The values of the other parameters 
are xo = 0, fco = 1.7A - , S = 1.1A, Ao = 16.1.APosition x and momentum k are measured in A and A -1 , respectively. Notice 
the strong suppression of interference at large values of momentum, both in (a) and (b). In case (a) only one of the two 
Gaussians interacts with the fluctuating phase shifter; moreover, the interference term in (E4j) depends on x and the oscillating 
part of the Wigner function is bent towards the negative a;-axis. 
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FIG. 4. Phase shift A in @ and distribution function w(A) in @, for different values of rf. (a) n = 1/2; (b) r 2 = 2/3; (c) 
r3 = 3/5; (d) r^ = 5/8; (e) r$ = 8/13; (f) foo = (V§ — l)/2. In each figure Ai = 2 (we set Ao = for clarity of presentation): 
above, phase shift A(t); below, distribution function w(A). Notice that, by increasing j (index of the Fibonacci sequence), 
the two frequencies become mutually "less rational," the phase shift A(t) becomes more irregular and its distribution function 
u>(A) more uniform. [The entropy behaves accordingly, increasing for j = 1, . . . , 5 (see Fig. BKa))]. Notice that the number of 
divergences of the distribution function increases with j; as shown in Appendix B, in the j = oo limit, the dist ribution becomes 
continuous with only one (logarithmic) divergence in A = and can be expressed as an elliptic integral (B.6). 
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FIG. 5. (a) Entropy ( p7| ) vs j (index in the Fibonacci sequence), (b) Decoherence parameter (Ft7j) vs j: case of a single 
Gaussian, (c) Decoherence parameter (|17|) vs j: case of a double Gaussian in a magnetic field. Notice that, while the entropy 
is an increasing function of j for j — 1, . . . , 5, the decoherence parameter displays a maximum at j — 3, both for a single and a 
double Gaussian. 



TABLE I. Entropy and Decoherence 



:i 


r l 


S 


e (single Gaussian) 


e (double Gaussian) 


l 


1/2 


1.6165 


0.52894 


0.59545 


2 


2/3 


1.7398 


0.53166 


0.62478 


3 


3/5 


1.7458 


0.53199 


0.63184 


4 


5/8 


1.9051 


0.53173 


0.62695 


5 


8/13 


1.9434 


0.53173 


0.62695 
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